function plot_m2d2(x,y)
clf

x0=0
x1=3

%a1 = x0:(x1-x0)/10:x1;


%for a1=.4
%g=1.5
lw=2;
fs=22;
d=x1/1000;

c1 = 2*(2-sqrt(2));
c2 = 2;

ylim([0,1])


plot([c1 c1], [0  1]);
hold on;
plot([c2 c2], [0  1]);
%UNTITLED2 Summary of this function goes here
%   Detailed explanation goes here
%plot(x,y,'LineWidth', lw);
%legend('','','f(0,0)', 'f(0,1)','f(1,0)', 'f(1,1)');
plot(x,y(:,1),'m',x,y(:,1),'ms','LineWidth', lw);
plot(x,y(:,2)+y(:,3),'r',x,y(:,2)+y(:,3),'ro','LineWidth', lw);
plot(x,y(:,4),'g',x,y(:,4),'gd','LineWidth', lw);

xlabel('\gamma', 'FontSize', fs);
ylabel('f', 'FontSize', fs);
set(gca,'FontSize',fs);

text(.84*x1, .06, 'f(0,1)+f(1,0)', 'FontSize', fs)  
text(.93*x1, .375, 'f(0,0)', 'FontSize', fs)  
text(.93*x1, .63, 'f(1,1)', 'FontSize', fs)  
%fs=16;
%text(c1/2 -.3, .95, 'B'' = \{(0,0) (1,0)\}', 'FontSize', fs)  
%text(c1 + (c2-c1)/2 -.41, .95, '\{(0,0) (1,0) (1,1)\}', 'FontSize', fs)  
%text(c2 + (x1-c2)/2 -.3, .95, 'B'' = \{(0,0) (1,1)\}', 'FontSize', fs)  


name = './eq';
print('-depsc',strcat(name, '.eps'));
%print('-dpdf',strcat(name, '.pdf'));
hold off;

end

